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The extraordinary radio galaxy MRC B1221— 423: probing 
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ABSTRACT 

We present optical spectra and high-resolution multi-wavelength radio observa- 
tions of the compact steep-spectrum radio source MRC B1221— 423 [z = 0.1706). 
MRC B1221— 423 is a very young (~ 10^ yr), powerful radio source which is un- 
dergoing a tidal interaction with a companion galaxy. We find strong evidence of 
interaction between the AGN and its environment. The radio morphology is highly 
distorted, showing a dramatic interaction between the radio jet and the host galaxy, 
with the jet being turned almost back on itself. H i observations show strong absorp- 
tion against the nucleus at an infall velocity of ^ 250 kms^^ compared to the stellar 
velocity, as well as a second, broader component which may represent gas falling into 
the nucleus. Optical spectra show that star formation is taking place across the whole 
system. Broad optical emission lines in the nucleus show evidence of outflow. Our ob- 
servations confirm that MRC B1221— 423 is a young radio source in a gas-rich nuclear 
environment, and that there was a time delay of a few xlOO Myr between the onset 
of star formation and the triggering of the AGN. 
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1 INTRODUCTION 

There is growing evidence of strong links between the trig- 
gering of powerful radio sources and gravitational interac- 
tions between galaxies. Many radio galaxies show signatures 
of tidal interact ions, such as tai l s, bri d ges, shells, and dou- 
ble n uclei (e.g. iHeckman et al.l Il98a : iBarnes fc Hernguistl 
11993 ). These tidal encounters can lead to large increases in 
the amount of material being fed to the central black hole, 
thereby triggering the strong radio emission. The same tidal 
event that dumps gas in the vicinity of the nuclear black 
hole can produce a burst of star formation, and hence a 
population of young stars. 

Further, the mechanical energy associated with rela- 
tivistic radio jets can have a profound effect on the formation 
and evolution of galaxies, by heating the interstellar gas to 
the point where stars can no longer form. This radio-mode 
feedback has been suggested as the mechanism which sup- 
presses star formation in massive galaxies; without such a 
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feedback mechanism, galaxy evoluti on models do not m atch 
the observed properties of galaxies IjCroton et al.ll2006l ). 

So both starbursts and the growth of the nuclear black 
hole are intimately associated with the merger history of 
the galaxy. However, key questions remain unanswered. If 
mergers are responsible for the formation of both AGN and 
starbursts, why do we sometimes see one without the other? 
When we do see both, we usually see a young radio source 
together with an older starburst, suggesting that there is a 
delay between the starburst and the onset of AGN activity. 
What is the reason for this delay? 

The compact steep-spectrum (CSS) sources are ideal 
laboratories for testing these connections. CSS sources are 
powerful but compact radio sources, with sizes 1-20 kpc, so 
that the radio source lies well within its host galaxy. Their 
radio spectra are steep, a < —0.5 {S„ oc v""), peaking below 
500 MHz. They have the same distribution in r edshift as the 
rest of the powerful radio source population fsee lO'Deall998l . 
for a review). The weight of evidence now is that CSS sources 
are small because they are young; if fuelling is sustained they 
are expected to evolve into classical large double-lobed radio 
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sources (|0'Dealll998l ). All CSS sources accessible to detailed 
optical investigations show evidence for recent interactions 
and/or mergers, and evidence for young stellar populations 
and galaxy- wide starbu rsts have been found in many of these 
objects (see lHoltl l2009. for a review). 

MRC B1221-423 is a member of a sample of 
southern radio sourc e s simi lar in power to 3C sources 
(|Burgess fc HunsteadI l2006al ). with a radio power of 
Pi.4GHz = 1.8 X 10^" WHz"^ The optical emission of the 
host galaxy is distorted, showing clea r signatures of inter - 
action, including tidal tails and shells l|Safouris et al.ll2003h . 
Multicolour images show knotty blue star- forming regions, 
while the likely merging galaxy to the south is also signifi- 
cantly bluer than the host galaxy (Johnston et al. 2005; see 
also F igure [T]). With a redshift oi z = 0.1706 (|Simpson et al.l 
ligoaj ). MRC B1221-423 is one of the nearest CSS sources; 
in the 3C sample, we need to go to z = 0.3 to find the first 
compact source as powerful as MRC B1221— 423. The radio 
source has double lobes which lie well within the envelope 
of the host galaxy (jSafouris et al.ll2003l ). The age of the ra- 
dio source, estimated from the power of the radio jets and 
total stored energy in the s ynchrotron plasma of the lobes 
IJRawlings k, SaundersHlGQll ). is only ~ 10^ yr. As a young, 
nearby, bright CSS source, MRC B1221-423 is therefore an 
ideal candidate for testing the relationship between the AGN 
and its environment. 

The goal of this study is to understand the relationship 
between the central AGN, which is the observational sign- 
post of a growing black hole, and the environment of the 
galaxy. MRC B1221— 423 provides a unique opportunity to 
measure both the age of the radio source and the star for- 
mation history of the galaxy, both of which are connected 
to the merger history of t he galaxy. 

In a previous paper (jJohnston et al.l [2005J, henceforth 
Paper I) we presented an analysis of multi-colour BVRIK 
images of MRC B1221— 423, which we used to analyse the 
change in stellar populations across the galaxy and its com- 
panion. We found evidence for three distinct stellar popu- 
lations, with ages of about 15 Gyr, 300 Myr and less than 
10 Myr. We concluded that the interaction with the com- 
panion galaxy triggered star formation in both galaxies, and 
after a substantial time delay, the infalling gas triggered ac- 
tivity in the central supermassive black hole. 

In order to better understand the connection between 
the radio source and its environment, we have undertaken 
further radio and optical observations of MRC B1221— 423, 
which we present in this paper. In sections [2] and [3] we 
describe the radio images and Hi absorption spectra of 
MRC B1221— 423, and in section |4] we describe the optical 
spectra. We have assumed a flat A-CDM cosmology with 
Ho = 71 kms"^ Mpc"\ ^m = 0.27 and f^A = 0.73 through- 
out. At the distance of the galaxy, 1 arcsec corresponds to 
a distance of 2.876 kpc. 



2 RADIO IMAGING 

We observed MRC B1221-423 at 1344, 2282 and 18 880.5 
-I- 19 648.5 MHz (20 cm, 13 cm and 12 mm wavebands) with 
the Australia Telescope Compact Array (ATCA), and also 
extracted unpublished, archival data at 1384 and 2496 MHz 
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Figure 1. V- band NTT image of the liost galaxy (from 
Ijohnston et aI.M 2005\ showing knotty features and tidal tails, 
strongly suggesting an ongoing interaction with the companion 
galaxy (arrowed). The 12 mm radio contours (Fig. [2]| are over- 
laid to show the scale of the radio source in relation to the galaxy. 
The two slit positions (Table[3]l arc shown; spectrum A was taken 
along the rotation axis of the galaxy (see § 14.11 1 , while spectrum 
B was taken along the line joining the galaxy to its companion, 
which is close to the axis of the radio source. 



from the ATCA online archivqj. Table [T] shows the details 
of these observations. 

Th e data w ere reduced using standard procedures in 
MIRIAD (jSault e t al. 1995). The effective bandwidth is either 
104 MHz (13 channels, each of width 8 MHz) or 56 MHz (14 
channels, each of width 4 MHz). However, at both 1384 and 
2282 MHz, not all of the band is usable because the ATCA 
suffers from self-generated interference at integer multiples 
of 128 MHz. Hence, there is a slight shift of the effective 
frequency in both cases (see Table [l|. 

PKS B1934— 638 was the primary calibrator at 20 and 
13 cm, except for the observations from 2003 November 
where it was not observed. In this case, PKS B0823— 500 
as well as the 2280 MHz data from 2003 August were used 
to determine the flux density scale; we estimate that the 
uncertainty in this procedure is a few per cent at most. At 
18.9 and 19.6 GHz, flux calibration was carried out using 
observations of Jupiter; we used the shortest baseline only, 
where Jupiter is not significantly resolved out. 

The deconvolution and imaging procedure at each fre- 
quency consisted of a number of iterations of CLEAN plus 
self-calibration. Only after several rounds of phase-only self- 
calibration was a full amplitude and phase solution at- 
tempted. To ensure that the source model derived from the 
CLEAN components was as simple as possible (especially in 
the 12-mm waveband), natural weighting was used while 
self-calibrating. When it was clear that there was no fur- 
ther improvement in the image quality, both naturally- and 

^ http://atoa.atnf.csiro.au 
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Figure 2. Total-intensity, uniformly-weighted contour map of 
MRC B1221-423 at 19.3 GHz, with the position of the optical nu- 
cleus shown in greyscale. The lowest contour is 0.75 mjy beam~^, 
with each successive contour corresponding to a doubling of the 
surface brightness. Vectors show the projected electric field at 
each pixel location (not corrected for Faraday rotation; pixel size 
0.1 arcsec X 0.1 arcsec); the scale bar in the bottom right-hand 
corner corresponds to 50 per cent fractional linear polarisation. A 
S/N cutoff of 5 in polarised intensity was used when plotting the 
polarisation vectors. The various components have been labelled. 
The synthesised beam is shown in the bottom right-hand corner. 



uniformly-weighted images were created from the final self- 
calibrated visibilities. Table [T] contains a summary of the 
angular resolution and RMS noise level in the vicinity of 
MRC B1221-423 in each total intensity image. The 18.9 
and 19.6 GHz images were restored using the same beam 
size (that of the 18.9 GHz data) to allow these images 
to be averaged together for further analysis (average fre- 
quency 19 264.5 MHz). We estimate that the astrometry in 
the millimetre- wavelength images is accurate to ~0.1 arcsec. 



2.1 Morphology 

At 8640 MHz, MRC B1221-423 is resolved into a north- 
south double structure with a projected separation of 4.3 
kpc (1.5 arcsec), as well as a third component to the east 
ISafouris et al.ll2003l ). However, the finer angular resolution 
at 19.3 GHz reveals that MRC B1221-423 has a particu- 
larly striking radio morphology akin to a string of pearls 
(Figure [2]). The southern jet has undergone a dramatic in- 
teraction with the host galaxy, being bent through a full 
180°. The projected length of the jet from the position of 
the optical nucleus to the peak of the furthest knot (knot 
F) is 8.6 kpc. On the other hand, the northern jet is much 
shorter: the projected separation is only 2.5 kpc. 



2.2 Spectral energy distribution 

Figure [3] shows the spectral energy distribution of 
MRC B1221-423; the fluxes are tabulated in Table [XT] Be- 
cause the source is only slightly extended at 20 and 13 cm, 
the 1344, 1382, 2280 and 2496 MHz flux densities were each 
determined by fitting an elliptical Gaussian with the miriad 



Figure 3. The observed-frame radio SED for MRC B1221-423, 
using all of the flux densities from Table |Xl] The 843, 4800 and 
4850 MHz data points are the non- weighted averages of the two 
individual measurements at each of these frequencies. A second- 
order polynomial has been fitted to the SED: ln(Si/) = 10.15 — 
0.01619 ln(!/) - 0.04295[ln(^)]2 with v in MHz and Sv in mJy. 
Though nearly all of the Parkes flux densities are affected by 
blending, the fit to the SED is not significantly altered if these 
data points are excluded from the fit (the two fits agree typically 
to within a few per cent). 



task IMFIT. However, it was found that the most accurate 
way to measure the flux density at millimetre wavelengths 
was to sum the flux in a polygon that enclosed the source. 
The uncertainty for each of the new ATCA measurements is 
dominated by the internal calibration uncertainty. While the 
20 and 13 cm flux densities are not affected by the weight- 
ing scheme used, we have quoted the 18.9 and 19.6 GHz 
flux densities from the naturally-weighted images, which are 
approximately three per cent higher than the values from 
the uniformly- weighted images. This will be due to some 
of the extended emission being resolved in the uniformly- 
weighted images because of the flner resolution. However, 
the naturally-weighted images themselves may be missing 
some extended flux, as, for example, can be seen by com- 
pa ring our results wi th th e ^^20 GHz flux densi ties reported 
bv lRicci et aD (J2006') and lMurphv et all (J2010l l. which were 
derived using lower-resolution ATCA data than ours. 

The observed-frame radio spectral energy distribution 
(SED) between 80 MHz and 22 GHz steepens with increas- 
ing frequency (Figure |3} . The spectral index a varies from 
-0.53 at 408 MHz to -0.88 at 22 GHz. Moreover, the 1.4 
GHz luminosity is 1.8 x 10^^ W Hz~^ for our assumed cos- 
mological model. Further low-frequency measurements are 
needed to determine accurately the turnover frequency. 



2.3 Linear polarisation 

ISafouris et al.l ([2003) investigated the polarimetric proper- 
ties of MRC B1221-423 at 4800 and 8640 MHz, finding 
that that the fractional linear polarisation m is as high 
as 35 per cent at 8640 MHz (resolution 1.4 arcsec x 0.63 
arcsec). We now extend this analysis with our new ATCA 
high-frequency data. In Figure [21 we have plotted polarisa- 
tion vectors at 19.3 GHz (resolution 0.57 arcsec x 0.34 arc- 
sec) . The vectors themselves were computed by first forming 
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Table 1. Log of ATCA radio observations. 







Observing 


Effective 


Integration 


Angular resolution" 


Beam 


RMS 


noise*^ 


Date 




frequency 


bandwidth 


time 


(arcsec X 


arcsec) 


PA*" 




(mjy beam^^) 


(UT) 


Array 


(MHz) 


(MHz) 


(h) 


N 


U 


N 


U 


N 


U 


2004 Jan 9 


6A 


1344 


104 


8.44 


14.05 X 9.58 


8.12 X 5.71 


-1?1 


-0?3 


0.45 


0.55 


1997 Jul 7 


6A 


1384 


104^* 


0.06 


13.48 X 10.02 


10.12 X 6.60 


-24?3 - 


-24? 8 


1.2 


1.4 


2003 Aug 23 
2003 Nov 22 


EW367 
1.5D 


2282 
2282 


56<^ 


3.37 
6.18 


14.19 X 13.64 


4.20 X 3.31 


16° 


16° 4 


0.42 


0.55 


1997 Jul 7 


6A 


2496 


104 


0.06 


7.45 X 5.60 


5.33 X 3.17 


-24?4 - 


-23? 9 


1.3 


1.8 


2003 Jul 27 


6D 


18 880.5 


104 


7.59 


1.09 X 0.67 


0.57 X 0.34 


9? 6 


8?0 


0.11 


0.18 


19 648.5 


104 


7.59 


1.09 X 0.67 


0.57 X 0.34 


9? 6 


8?0 


0.11 


0.18 



■^N — natural weighting; U - uniform weighting. 

'Measured north through east 

'^RMS noise level in the vicinity of MRC B1221— 423 in the total intensity image 

■^One of the 8 MHz channels is flagged due to self-interference at 1408 MHz. The effective frequency is shifted slightly to 1382 MHz. 

"^One of the 4 MHz channels is flagged due to self-interference at 2304 MHz. The effective frequency is shifted slightly to 2280 MHz. 



Table 2. 19.3 GHz flux density, fractional polarisation and po- 
larisation position angle measurements in the uniformly-weighted 
image (Fig. [2)l. The measurements were made at the position of 
the peak pixel in Stokes / for each of the source components. 



Component 


Si 


m 


X 




(mJy beam~^) 


(%) 


(°) 


A 


89.9 


13.0 ±0.2 


52.1 ±0.3 


B 


1.5 


< 20° 




C 


48.4 


1.9 ±0.2 


59 ±3 


D 


42.0 


5.7 ±0.3 


14.3 ±1.2 


E 


11.8 


14.7 ±0.9 


-67.0 ±1.6 


F 


11.9 


14.7 ±0.9 


-1.3 ±1.6 



°3(T upper limit. 



Stokes Q and U images at both 18.9 and 19.6 GHz from the 
self-calibrated visibilities; these images were restored with 
the same synthesised beam size. Averaged images were then 
produced for Stokes Q and U at 19.3 GH z, from which a to- 
tal polarised intensity (P — ^JQ'^ +U^) map, corrected for 
Ricean bias, was also formed. The average RMS noise level 
in the 19.3 GHz Stokes Q and U images (ctqu) in the vicin- 
ity of MRC B1221-423 is about 0.1 mJy beam"\ Note that 
bandwidth depolarisation is negligible for MRC B1221— 423 
in the 18.9 and 19.6 GHz bands. 

High fractional linear polarisation is also observed at 
millimetre wavelengths, particularly at the eastern edge of 
the source, where the mean fractional polarisation is ~50 per 
cent, with some pixels having m > 60 per cent. The fractional 
polarisation is also significant at the southern edge of the 
source (30-40 per cent), as well as between knots D and 
E (~25 per cent on average). Over the entire source, the 
median fractional polarisation of the pixels is about 14 per 
cent. Further polarisation statistics are shown in Tabled 

The orientation of the projected magnetic field with re- 
spect to the jet path varies significantly across the source. 
Knots A and C have similar polarisation position angles even 
though there is a significant difference in the degree of po- 
larisation; in these regions the magnetic field is neither par- 
allel nor perpendicular to the jet path. While the magnetic 
field becomes increasingly perpendicular to the jet path at 
knot D, it then becomes strikingly parallel to the jet between 
knots D and F where the extreme bending occurs. The mag- 
netic field is also oriented circumferentially along the outer 



contours in this latter region, which was also observed i n 
the lower-resolution 8640 MHz data dSafouris et al.ll2003l ). 
In complete contrast, knot F exhibits magnetic fields that 
are closely perpendicular to the jet path. 

Note that after correcting for Faraday rotation (see be- 
low), the intrinsic and observed polarisation position angles 
are estimated to agree to within < 5° except for knot A, 
where the vectors need to be rotated anti-clockwise by ~10- 
15°. 



2.4 Rotation measure synthesis 

Faraday rotation measures (RMs) provide an important in- 
sight into the properties of the magnetoionic environment 
along the line of sight towards a radio source. Tradition- 
ally, the RM is determined from a linear fit between the 
polarisation position angle x s-nd A^. However, this method 
suffers from several problems, such as position angle ambi- 
guities. A superior Fourier-based technique, known as Fara- 
day rotation measure synthesis, overcomes the problems 
associated with x~'^ fit ting. Det ailed descriptions of the 
method can be fo und in I Brent i ens & dc Bruyn (2005) and 
iHeald et al.l (|2009|). In brief, the complex linear polarisation 
P(A^) — Q(A^) + iU{\^) and complex Faraday dispersion 
function F(<?!)) are related via a relationship that is very sim- 
ilar to a Fourier transform. The Faraday depth (j), which is 
a more general quantity than the RM, is defined as 



= 0.812 



ricB ■ dl 



(1) 



where L is the path length through the magnetised plasma 
in pc, rie is the thermal electron density in cm~^, and B is 
the magnetic field in fiG. Measurements of P in a number 
of channels across one or more bandpasses can thus be used 
to form an 'RM data cube'; a slice through this cube at a 
particular spatial position gives a Faraday depth spectrum 
(|P(</>)1 versus (j)). 



2.4-1 Low-frequency RM synthesis 

First, we performed RM synthesis at 20 and 13 cm using 
three different approaches: 1344 MHz alone (13 x 8 MHz 
channels), 2280 MHz alone (13 x 4 MHz channels) and by 
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combining the 1344 and 2280 MHz bands. There is no sig- 
nificant advantage in also including the 1382 and 2496 MHz 
observations because aqu is much higher than at 1344 and 
2280 MHz (~0.12 and ~0.20 mjy beam'^ respectively, 
in the two latter cases). Each Stokes Q and U individual 
channel image was restored at the resolution of the 1344 
MHz data (8.12 arcsec x 5.71 arcsec; PA = -0?3). Because 
MRC B1221— 423 is only slightly extended at these frequen- 
cies, we formed Faraday depth spectra at the position of 
peak intensity in Stokes I. Though we cannot spatially re- 
solve the individual components of the source, RM synthe- 
sis is still useful: the A^ coverage is greater than at higher 
frequencies, which reduces the uncertainty in the Faraday 
depth because the resolution in ^-space is much higher. 
Therefore multiple Faraday depth components can poten- 
tially be resolved in 0-space even though we cannot spatially 
distinguish the different Faraday screens. 

Inverse- variance weighting was applied to each channel 
when forming the Faraday depth sp ectra. A spectral index 
corre ction was also implemented (see lBrentiens fc de Bruvnl 
|2005|). The spectra were CLEANed to remove the sidelobes 
that are present due to incomplete A^ coverage; the response 
function is referred to as the rotation measure spread func- 
tion (RMSF). A Gaussian 'restoring bea m' was used when 
constructing each CLEANed spectrum fsee lHeald et al.ll2009l . 
for further details); the resolutions in i/i-space are 591, 5123 
and 74 rad m"^ FWHM at 1344 MHz alone, 2280 MHz 
alone, and 1344 -f 2280 MHz combined, respectively. 

Both the 1344 and 2280 MHz spectra are found to con- 
sist only of a single point-like Faraday depth component. 
However, the component shifts position between the two 
frequencies: the position of the peak is —35 ± 4rad m~^ 
at 1344 MHz, but -220 ± 30 rad m"^ at 2280 MHz. The 
uncertainties take into account both the noise statistics and 
the estimated residual instrumental polarisation. The com- 
ponents are bright enough such that the Ricean bias is not 
significant. 

The shift in the position of the peak suggests that there 
is an unresolved distribution of components whose relative 
strengths vary with frequency. Indeed, multiple components 
are found when the 1344 and 2280 MHz data are combined 
to form a higher-resolution spectrum (Figure Q. Apart from 
the bright peak at —40 ± 2rad m~^, there are weaker com- 
ponents above lOcr (about 1.3 mJy beam~^ RMSF~^) at 
-1-270 and —510 rad vn~^ , the latter being extended over 1.6 
beamdwidths. However, these results must be interpreted 
with caution: the gap between the 20 and 13 cm bands in- 
troduces large sidelobes into the RMSF, which complicates 
the deconvolution process (though see the results at higher 
frequency below). 

The Faraday depth of the brightest component in Fig- 
ure |4] is ve ry similar to th e mean RM of —35 rad m~^ cal- 
culated by ISafouris et al.l (|2003) using the traditional lin- 
ear fitting_a££roach_at 4800 -I- 8640 MHz. As discussed 
bv ISafouris et al.l ([2003), this result is consistent with the 
Galactic foreground RM in this region of the sky. We were 
also able to construct Faraday depth spectra at 1344 MHz 
for two nearby sources in the field: SUMSS J122328-423436 
and SUMSS J122316-423922. Again, the derived Faraday 
depths are consistent with a Galactic magnetoionic screen 
being a significant, if not the main source of Faraday rota- 
tion. 




Figure 4. CLEANed 1344 -|- 2280 MHz Faraday depth spectrum 
for MRC B1221— 423 at the position of peak intensity in Stokes 
/. The dotted line indicates a Faraday depth of zero. The Ricean 
bias has not been removed. The most prominent peak can be 
associated with the foreground Galactic RM. 



2.4.2 High-frequency RM synthesis 

We also attempted RM synthesis at higher frequencies by 
combining re-analysed 8640 MHz data (uniform weighting; 
resolution 1.14 arcsec x 0.72 arcsec; PA = —9.4°) with our 
millimetre-wavelength observations (naturally weighted but 
restored at the resolution of the 8640 MHz data) ; there are 
a total of 39 channels, each of width 8 MHz. Though the 
angular resolution is superior at these frequencies, the A^ 
coverage is poorer: the FWHM in (/)-space is about 2600 rad 
m~^ for inverse-variance weighting, and, like at 1344 -I- 2280 
MHz, the RMSF contains large sidelobes because of the gaps 
between the bands. Because it is very easy for the peak of 
i^^((?!))| to be shifted to the position of a nearby sidelobe for 
the lower S/N pixels, we only measured Faraday depths for 
the pixels with peak polarised flux densities greater than 
1 mJy beam~^ RMSF~^, after applying a spectral index 
correction as above. 

Figure [5] shows a 3 cm -I- 12 mm Faraday depth map for 
MRC B1221— 423. Each pixel has only one component 
(though the sidelobes significantly complicate the search 
for secondary components), and none of these components 
are extended. The Faraday depth varies from —1260 to 
-1-330 rad m~^, with a steep gradient from north to south: 
the mean Faraday depth for the top half of the source is 
^ —800 rad m~^, while the mean for the bottom half is 
~ —30 rad ra~^ . The latter value is close to the Galactic 
foreground RM (see above) . The steep gradient in strongly 
suggests that the very large Faraday depths associated with 
the northern lobe are a consequence of a magnetoionic screen 
that is local to the source. The rest-frame Faraday depth will 
be a factor of (l-l-z)'^ = 1.37 times higher than the observed- 
frame value after the Galactic contribution is removed from 
the latter. Assuming that the foreground Faraday depth is 
~ —40 rad m~^, then the rest-frame Faraday depth at the 
position of peak intensity in Stokes / for the northern lobe is 
about —1100 rad m~'^. From Equation[T] large values of phi 
are a consequence of one or more of the following along the 
line of sight: (1) a dense medium, (2) a highly magnetised 
medium, and (3) a long path length. The radio observations 
alone are not sufficient to disentangle the relative contribu- 
tions. 

To investigate further the magnetoionic environment of 
MRC B1221— 423, we have also constructed a map of the de- 
polarisation ratio »ni9264.5/?728640 (Figure [6ll. The northern 
lobe is significantly depolarised with decreasing frequency, 
implying small-scale inhomogeneities in the Faraday screen. 
This is likely to explain why there is not a prominent 4> 
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8.6 + 18.9 + 19.6 GHz 



-42:35:34.29 




Figure 5. 8.6 + 18.9 + 19.6 GHz observed-frame Faraday depth 
map for MRC B1221— 423 (angular resolution resolution 1.14 arc- 
sec X 0.72 arcsec; beam PA —9.4°; pixel size 0.23 arcsec X 0.23 
arcsec). The total-intensity radio contours are at 8640 MHz; the 
lowest contour is 2 mjy beam~^ and each successive contour cor- 
responds to a doubling of the surface brightness. The colour bar 
has units of radm~^. The typical uncertainty in the Faraday 
depth for a given pixel is about 30 rad m~^. The letters label the 
positions of the components seen in the 19.3 GHz map (Fig.[2]l. 




I 






Figure 6. Depolarisation ratio map between 8.6 and 19.3 GHz. 
The radio contour scheme is the same as in Figure [5] A S/N 
cutoff of 5 in polarised intensity was used at both 8.6 and 19.3 
GHz when constructing the map. 



3 RADIO SPECTROSCOPY 

The Hi observations of MRC B1221-423 were done with 
the ATCA for 2 x 12 hrs on 2004 January 9 and 2004 July 
17 using the 6-km array-configuration. All observations were 
done with 16 MHz bandwidth and a central frequency (for 
the first IF) set to 1214 MHz. The first observation used 256 
channels and two IFs while the second observation used 512 
channels and only one IF. 

We used PKS B1934— 638 as bandpass and flux cal- 
ibrator. To monitor the system gain and phase changes, 
PKS B1215— 457 was observed as a secondary calibrator for 
ten minutes every hour. 

For the data re duction and visualisation we used the 
MIRIAD and KARMA (|Goochlll996r ) data reduction packages. 
For the line observations, after flagging and calibration, we 
separated the continuum from the line data in each indi- 
vidual data set by fitting either a first or a second order 
polynomial to the line-free channels. 

The continuum was subtracted using uvlin by making 
a second-order fit to the line-free channels of each visibil- 
ity record and subtracting this fit from the spectrum. After 
cleaning, Hanning smoothing was applied to each line cube. 
The final cube was derived by combining the datasets for the 
two observations. The beam size in the resulting uniformly 
weighted cube is 9'.'6 x 6'.'4, PA = —1°. The rms noise in 
the Hanning cube is 0.77 mJy beam~^ per channel. Velocity 
resolution is 26 kms~^ after Hanning. 

H I in absorption was clearly detected. The absorption, 
shown in Figure [T] consists of a deep and relatively nar- 
row component with a broad wing. The narrow compo- 
nent (FWHM ~ 90 kms"^) has a peak of -40.5 mJy cen- 
tred on ~ 51370 kms^^. The optical depth of this com- 
ponent is r ~ 0.02. Compared to the systemic velocity de- 
rived from the stars (51118 ± 15 kms~^; see ij |4.2|l . this ab- 
sorption component is redshifted by ^ 250 kms~^. The 
broad wing (estimated full-width at zero intensity (FWZI) 
~ 400 kms~^) is blueshifted with respect to the narrow com- 
ponent by ~ 150 kms~^. It covers a velocity range extend- 
ing to ~ 51100 kms~^, therefore reaching velocities slightly 
blueshifted compared to the systemic velocity of the galaxy. 
The peak absorption of this wing has an optical depth of 
T ^ 0.003 (see Figure [T]!. 

No blueshifted component (which could be connected 
with outfiowing gas) was detected, to an optical depth limit 
of r ~ 0.001 (3cr). This is comparable with the optical depth 
detected in broad absorption co mponents associated with 
outflows in other radio galaxies (JMorganti et al.l 120051 ). so 
we would have been able to detect such outflows if they 
were present. 



4 OPTICAL OBSERVATIONS 



component at ^ —1000 rad m~^ in Figure |4] On the other 
hand, the southern jet is not significantly depolarised be- 
tween 19264.5 and 8640 MHz; the fractional polarisation in 
fact increases with decreasing frequency in some areas. 

It is possible that we are seeing an orientation effect, 
with the radio source fore-shortened and the northern lobe 
seen through the galaxy. We discuss this possibility in Sec- 
tion [5] 



4.1 Observations and data reduction 

Optical spectra were obtained on 2004 May 14-16 using the 
Double Beam Spectrograph on the ANU 2.3-m telescope at 
Siding Spring Observatory (Table[3ll . The detectors were two 
SITe 1752 x 532 CCDs with 15-^m pixels. For the first night, 
a dichroic filter with a cross-over wavelength of 6300 A was 
used to split the light into the two arms of the spectrograph. 
Two gratings with 600 grooves mm~^ were used, giving a 
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Table 3. Journal of optical spectroscopic observations of 
MRC B1221-423. Two different set-ups were used, on the first 
night and remaining nights respectively. The columns show the 
date of the observation, the wavelength range and the resolution, 
the exposure time, the position angle of the slit on the sky, and 
the median seeing of the observations. The final column shows 
the slit position (Fig. [ij. 



51000 51500 52000 

Figure 7. Hi absorption against the nucleus of MRC B1221— 423, 
showing large amounts of redshifted (infalling) material compared 
to the systemic velocity (dashed line) . There is also a broad wing 
of absorption extending blueward of the narrow absorption com- 
ponent. The thick red line shows the sum of two gaussians fitted 
to the profile; the individual components are shown as the dot- 
ted blue lines. The vertical dashed line with error bar shows the 
stellar velocity {cz = 51118 ± 15 kms~^, see ij 14. 2. 31 1. while the 
vertical dash-dotted line shows the mean omission-line velocity 
{cz = 51166 ± 30 kms-l, see S l43)l . 



wavelength resolution of 2.2 A over the wavelength range 
420 0-6120 A in the blue and 6270-7130 A in the red. 

ISafouris et al.l (120031 ) used spectra taken at three dif- 
ferent position angles to determine the angle of maximum 
gas rotation velocity; the (projected) axis of rotation was 
assumed to be perpendicular to this. The slit was oriented 
along this rotation axis, at a position angle of 118° (slit posi- 
tion A; Figure[l|. This set-up gave us high-resolution spectra 
with minimum rotational broadening of the line profiles. 

For the remaining two nights, a plane mirror was used 
in place of the dichroic, so all hght was sent to the blue arm 
of the spectrograph. A 300 grooves mm~^ grating gave us a 
wavelength resolution of 4.3 A over the range 3700-7500 A. 
The slit was oriented at a position angle of —18° so as to 
place both the galaxy and its interacting companion on the 
slit (slit position B; Fig.[l]). This set-up gave us slightly lower 
wavelength resolution but continuous wavelength coverage 
of the blue region (rest-frame 3160-6400 A), which is useful 
for age-dating stellar populations. 

The spatial scale on the detector was O'.'91/pixel, corre- 
sponding to a linear scale of 2.6 kpc per pixel. On all three 
nights, a slit-width of l'.'5 was used, and pairs of 1800 s ex- 
posures were bracketed with NeAr a rc-lamp exposures. 

The IRAF software suite l|Todvl|l986 ) was used to re- 
move the bias and pixel-to-pixel gain variations from each 
frame. Cosmic ra y events wer e removed using the tech- 
nique described b y lCrokd (|l995l ). as implemented in FIGARO 
l|Shortridgell 19 93). The spectra were straightened in FIGARO 
so that the dispersion ran exactly along rows of the image, 
then a two-dimensional wavelength fit was performed to the 
arc images by fitting a third-order polynomial to the arc 



UT date 


A range 


FWHM 


Ccxp 


PA 


Seeing 






(A) 


(A) 


(ks) 


(°) 


(arcsec) 




2004 May 14 


4200-6120 


2.2 


14.4 


118 


1.4 


A 




6450-8000 


2.2 


12.6 


118 


1.4 


2004 May 15 


3700-7500 


4.3 


23.4 


-18 


1.6 


B 


2004 May 16 


3700-7500 


4.3 


23.4 


-18 


1.5 



wavelengths as a function of pixel number, for each row of 
the image. These wavelength solutions were copied to the 
object images, interpolating between the bracketing arc ex- 
posures, and the data were rebinned so the wavelength-pixel 
relation was linear and uniform across the image. The sky 
background was subtracted from each image, by making a 
linear fit to the sky on either side of the galaxy, choosing 
regions well outside the wings of the galaxy profile. The 
spectra were corrected for atmospheric extinction, the tel- 
luric absorption features were removed by comparing with 
the spectrum of a smooth-spectrum standard taken at sim- 
ilar airmass, and the spectra were flux-calibrated. We con- 
structed two grand-sum spectra of the nucleus - one for 
each set-up (Table [31 - by aligning the individual spectra 
and summing, then correcting the su mmed spectrum for t he 
Galactic reddening Av = 0.331 mag (|Schlegel et al.lll998l ). 

These spectra are shown in Figure [S] They show nar- 
row emission lines from ionised gas, on top of a continuum 
dominated by hght from stars. In the following sections, we 
analyse these two components separately. 



4.2 Stellar continuum 

4.2.1 D4000 strength 

The continua of many radio galaxies show UV excesses com- 
pared with normal elliptical galaxies. This excess may arise 
either from scattered nuclear hght or from a population 
of young stars. In Figure |9l we show the UV excess in 
MRC B1221-423, as measured by the stren gth of the 4000 A 
break We use D'(4000) as defined by iTadhunter et all 
(2002), which is the ratio of the fiux between 4150-4250 A 
(above the 4000 A break) to the fiux between 3750-3850 A 
(below the break). T his is a version of the original index 
D(4000) proposed by iBruzual (|l983l ). modified to exclude 
emission lines which are often present in the spectra of radio 
galaxies. Old stellar populations have large values of this in- 
dex, while smaller values indicate a less pronounced 4000 A 
break, and hence the presence of more UV light than ex- 
pected from a passively evolving stellar population. 

Figure [9] shows that there is an excess of UV fiux in all 
regions of MRC B1221— 423 compared to a typical ellipti- 
cal galaxy. Nowhere does the observed D'(4000) approach 
the value of 2.3 t ypical for a passively- evolving galaxy of 
solar metallicity ( Tadhunter et al.l 12005^ ; the mean value 
over the ga laxy (from the regions labelled Gl to G2) is 
1.48 ± 0.05. ITadhunter et all l|2002l ) measured D'(4000) for 
a sample of radio galaxies, and found values in the range 
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Figure 8. Spectrum of the nucleus of MRC B1221— 423, extracted from the central 2''7, with emission features identified. The top 
spectrum (spectrum A) is from the 2004 May 14 blue and red data, with the slit along the rotation axis of the galaxy; the bottom 
spectrum (spectrum B) is extracted from the combined data from 2004 May 15—16, with the slit along the line joining the galaxy to its 
companion. Spectrum A is offset in the y-direction by 4 X 10~^® ergs~^ cm~-^ A for clarity. 




Figure 9. Spatial profiles of the fiux along the slit in spectrum A 
in the rest wavelength ranges 3750—3850 A (cross-hatching) and 
4150-4250 A (hatching). The ratio of the two, D(4000), is shown 
as the solid line (right-hand scale); the dashed horizontal line 
shows D(4000) = 2.3, which is the value expected for a 12.5 Gyr 
elliptical galaxy of solar metallicity. Sub-solar metallicity {Z = 
0.4 Zq) lowers this ratio to D(4000) = 1.8. The ratio measured 
for MRC B1221— 423 is consistently below this, indicating the 
presence of a blue spectral component. The horizontal bars at 
the top of the plot indicated the regions extracted for the nuclear 
spectrum N, the companion galaxy C, and two regions away from 
the nucleus Gl and G2. 



~ 1.2-2.6. iTadhunter et al.l (|2005l ) found clear UV excesses 
in a sample of three powerful radio galaxies, and interpreted 
tiiis as evidence of star formation associated with the radio 
source. The fact that we see a similar UV excess throughout 
MRC B1221— 423 suggests that the star formation associ- 
ated with the radio source is galaxy-wide. 



4.2.2 Modelling the continuum 

Next, we model the shape of the continuum spectral en- 
ergy distribution in detail. The spectrum from 2004 May 
14 (spectrum A) had different resolution, wavelength cover- 
age, and position angle to the grand-sum spectrum from the 
other two nights (spectrum B; Table |31), so we performed 
the following analysis for both summed spectra separately. 
We fitted three separate models to the continuum: 

(i) a single-age stellar population, allowing the age to vary 
(ii) an old stellar population plus a power-law continuum 
representing the AGN, with F\ on X"' , allowing a to vary 

(iii) an old stellar population plus a single-age young stel- 
lar population, allowing the age of the young stars to vary. 

We used published iso chrone spectral synthesis models 
from the GISSEL96 library l|Bruzual fc Charlojl 19931 ) for the 
stellar populations. Both the "old" and "young" components 
were single-age stellar populations, representing instanta- 
neous bursts of star formation. We used a Salpeter initial 
mass function with m ass limits of 0.1 and 100 Mq and solar 
metallicity, using the iGunn fc Strvkeij (119831 ) stellar spec- 
tral atlas. For the old population, we used an instantaneous 
burst model with an age of 15 Gyr; the exact age of the 
"old" population is not critical, since the spectrum changes 
little at these ages. For the young stellar population, we 
used instantaneous burst models with ages 10, 15, 20, 50, 
80, 100, 150, . . . 5000 Myr, giving us 14 models in total. We 
assembled the model spectra from the GISSEL96 library, and 
binned them to the same resolution as the observed spectra. 

The spectra (observed and model) were binned into 20 A 
bins, after which we removed bins corresponding to emission 
lines, residual night-sky lines, cosmic rays and atmospheric 
absorption bands. This resulted in 83 wavelength bins for 
spectrum A and 113 bins for spectrum B. We measure the 
flux in these bins for both our observed spectrum and mod- 
els. The models were matched to the observed flux in a bin at 
(rest-frame) 4165 A; the normalisation was allowed to vary 
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Table 4. Fits to the stellar continuum, in the nucleus and other regions of the galaxy. Fits are shown to spectrum A, taken along the 
rotation axis of the galaxy, and spectrum B, along the line joining the galaxy with its companion (lower resolution; TablefSJI. The three 
different models fitted to each spectrum are shown in column 2; column 3 shows the percentage of light at 4156 A contributed by the blue 
component (where present); column 4 shows either the index of the power-law component or the age of the single-age stellar population 
component; column 5 shows the reduced x^ of the best fit. Ranges on the parameters indicate l-cr errors on the best fit. The lower half 
of the table show the results of fits to the companion galaxy and to the main galaxy outside the nucleus. 



Spectrum 




Model 


% blue component 


Parameter 


xi 


Nucleus (A) 




Single age 


- 


2.2-2.7 Gyr 


4.18 






15 Gyr -|- power-law 


44 ±5 


a = -0.54 ±0.55 


2.14 






15 Gyr + young stars 


42 ±1 


< 23 Myr 


1.35 


Nucleus (B) 




Single age 


- 


1.6-3 Gyr 


2.10 






15 Gyr + power-law 


43 ±5 


a = 0.14 ±0.35 


3.51 






15 Gyr + young stars 


43 ±1 


75-105 Myr 


1.63 


Companion (B) 


15 Gyr + young stars 


81 ±3 


150-200 Myr 


0.86 


Main galaxy (A) 


5; 5; 


64 ±4 


< 110 Myr 


0.79 


Main galaxy 


(B) 


" " 


64 ±4 


140-180 Myr 


1.10 



between 75 per cent and 125 per cent of the flux of the ob- 
served spectrum, to allow for possible differences in slope 
between the model and the observed spectrum. Varying the 
fraction of the blue component, the normalisation, and the 
power-law index or age of the young stars, the best-fitting 
model was found by calculating the x^ difference between 
the observed spectrum and each model spectrum. We as- 
sumed errors of ±4 per cent in each wavelength bin for the 
X^-fitting. 

We began by fitting to the spectrum of the nucleus of 
the galaxy, extracting the central three rows (region N in 
Figure O, which correspond to the central 2.7 arcsec, or 
7.8 kpc. The rest-frame wavelength range over which the 
fitting was done was 3600-6100 A for spectrum A, and 3300- 
6100 A for spectrum B. The results of this fitting are shown 
in Table |4] No single- age spectrum described the continuum 
well; the nominal best-fit age of 2 Gyr systematically under- 
estimates the flux in the red. Spectrum A was reasonably 
well modelled by the addition of a power- law with a = — 1, 
though the same model did not fit spectrum B well. The ad- 
dition of a population of young stars provided a much better 
fit to the spectrum. Both spectra A and B are well fitted by 
a model where ~ 43 per cent of the light at (rest-frame) 
4200 A is contributed by a young stellar population with an 
age of < 100 Myr. Spectrum A requires that the young pop- 
ulation have an age of less than 23 Myr, while spectrum B 
was adequately fitted with models having ages between 75 
and 105 Myr. The fit to spectrum A is shown in Figure [TOl 

Thus the continuum in the nuclear region is not well 
described by either a single-age population, or a combina- 
tion of old stars plus a power-law component. Instead, we 
require at least two stellar populations with different ages 
to reproduce the spectrum adequately. 

We repeated this analysis with spectra extracted from 
other positions along the slit. We extracted the spectrum 
at the position of the companion galaxy; we also extracted 
regions of the main galaxy away from the nucleus (Fig. llOf) . 
The spectrum of the companion galaxy was extracted from 
spectrum B, including five rows centred 9 arcsec south from 
the nucleus (region C in Fig. |9]). For slit position B, we 
could extract spectra from the southern and northern por- 
tions of the main galaxy separately. We extracted three rows 
(2.7 arcsec) beginning ±2.25 arcsec either side of the nucleus 
(regions Gl and G2 in Figure |9} and fitted the two spectra 



separately; the results for both northern and southern spec- 
tra were consistent. For slit position A, the signal-to-noise 
ratio of spectra extracted separately was insufficient, so we 
extracted the same regions but added the spectra together 
before fitting. As can be seen in Figure [TOl the spectrum of 
the galaxy outside the nucleus shows a much younger stel- 
lar population than the nucleus, with much stronger Balmer 
lines and a smaller 4000 A break. 

To quantify this, we again fitted the observed spectra 
with two single-age populations. The results of this fitting 
are shown in Table]?] The spectrum of the companion galaxy 
is well-described by a model where 80 per cent of the light is 
contributed by a 200 Myr old stellar population. The main 
galaxy away from the nucleus is best fitted with a slightly 
younger second population, with ages ~ 100 Myr. No sig- 
nificant difference was found between spectra obtained on 
either side of the nucleus. 

We investigated including internal reddening as an- 
other free parameter; the fits to the nuclear spectra were 
marginally improved by the inclusion of 0.4 mag of extinc- 
tion, but this did not significantly change the fitted ages of 
the underlying populations. The fit to the companion galaxy 
also indicates a reddening of Av ~ 0.4 mag, but the main 
galaxy outside the nucleus did not demand any reddening. 
These reddenings are significantly lower than the values we 
derived in Paper I, where we modelled the pixel colours us- 
ing single-age populations. 

This is broadly consistent with the results of Paper I, 
based on modelling of pixel colours, where we found an old 
population in the 'outskirts of the galaxy, an intermediate- 
age population (r ~ 300 Myr) in the companion and tidal 
tail, and a young population (r < 10 Myr) in the nucleus 
and blue "knots". Modelling the spectra enables us to best 
disentangle the contribution of the different populations, al- 
though our spatial resolution is not as good as in the images. 
Our fits to the pixel colours in the nucleus in Paper I were 
quite poor, indicating that our assumption of a single-age 
population was not a good one. The two-age fit to the spec- 
trum in this work clearly gives a better description. 

Thus we have found further evidence of starburst activ- 
ity throughout the whole MRC B1221-423 system, with 
stars of different ages present in different regions. The 
youngest stars are in the nucleus, with slightly older popula- 
tions in the companion galaxy and the main galaxy outside 



© 0000 RAS, MNRAS 000, 000-000 



10 H. M. Johnston et al. 




4000 5000 6000 4000 5000 6000 4000 5000 6000 
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Figure 10. The best two-component model fits to the continuum of the nucleus, the galaxy outside the nucleus (Gl, north of the 
nucleus), and the companion galaxy; the emission lines were not included in the fits. The apertures used to extract the spectra are shown 
in Figure |9] Note that the j/-scale of the left panel is double that of the other two panels (right-hand scale). The parameters of the fits 
are shown in Table |4] The model spectrum (thick black line) is shown on top of the observed spectrum, together with the individual 
components, the old population (red) and the young population (blue). 



of the nucleus. The stellar models enable us to determine 
what fraction of the mass of the galaxy comes from each 
population. The mass fraction in young stars is not large: 
the young population in the nucleus with age < 75 Myr, 
which contributes nearly half the light at 4200 A, contains 
less than 1 per cent of the mass. The highest mass fraction 
is found in the companion galaxy, where the 170 Myr pop- 
ulation contributing 80 per cent of the light makes up 5 per 
cent of the mass of the galaxy. 



i.2.3 Velocities 

We determined the rotation curve of the stellar popula- 
tion by cross-correlating the spectrum extracted at differ- 
ent positions along the slit with a synthetic spectrum from 
ICoelho et al.l (2005). We extracted non-overlapping 2-pixel 
spectra along the slit in spectrum A, the spectrum with 
the highest resolution, and cross-correlated with a synthetic 
spectrum of an F5V star (with Teff = 6750 K, log*; = 4.5), 
chosen as the best match to the continuum away from the 
nucleus of the galaxy. We cross-correlated the region from 
rest-frame 3845-4800 A, chosen to avoid the bright emis- 
sion lines, using the FIGARO rou tine scross, which uses the 
method of lTonrv fc Davia (|l979 ). The resulting stellar rota- 
tion curve is shown as the solid line in Figure [TS] 

The redshift of the nucleus, which we take as the "true" 
redshift of the galaxy, is 0.17053 ± 0.00005 {cz = 51118 ± 
15kms~^). 



4.3 Emission line gas 

The spectrum of MRC B1221— 423 shows numerous strong 
emission lines, characteristic of powerful radio sources. The 
emission lines enable us to probe the physical conditions 
and kinematics of the gas in the galaxy. Previous studies of 
CSS sources have found widespread evidence of oi itflows and 
disturbe d kinematics in these systems (see e.g. iHolt et al.l 
[2003: Emonts et al.ll200^ . 

We subtracted the best fit to the continuum from each 
of Spectrum A and spectrum B in order to model the emis- 



Table 5. Emission line measurements, from the rest-frame spec- 
trum B with the model stellar continuum subtracted. 



Line 



(xlQ- 



Flux 
erg s~ 



FWHM 

(km s"-*-) 



[On] 3727 


37.0 ± 1.5 


660 ± 40 


[Nein] 3869 


5.0 ±0.6 


730 ± 70 


Hen 4686 


< 0.12 


- 


H/3 4861 


6.1 ±0.4 


470 ± 40 


[O ni] 4959 


8.6 ± 0.2 


600 ± 30 


[Oni] 5007 


22.2 ± 0.2 


600 ± 30 


[O i] 6300 


7.0 ± 0.4 


680 ± 30 



sion lines. We used the fit comprising two single-age popu- 
lations, as determined in ij |4.2.2l and Table |4] 

Table [5] shows the measured parameters of the emission 
lines. The line centre and width (FWHM) were determined 
by fitting a Gaussian to the subtracted spectrum; the line 
flux was obtained by direct summation of the flux over a 
wavelength range four times the FWHM, except where such 
a width would also include neighbouring lines. 

A combined fit to the emission lines of the blue Spec- 
trum A gives a heliocentric redshift of 2: = 0.17068 ±0.00010 
{cz — 51168 ± 30kms~^); the red spectrum gives 2 = 
0.17087 ± 0.00013 {cz = 51225 ± 40kms-^). The emis- 
sion lines in Spectrum B give a redshift consistent with 
those in blue Spectrum A: z = 0.17067 ± 0.00010 {cz = 
51166 ± 30 kms"^). 

4-3.1 Kinematics of the emission lines 

CSS sources often show broad, asymmetric emission lines 
in the nuclear region. We extracted the spectrum at the 
position of the nucleus by selecting the central two rows 
from each two-dimensional spectrum, corresponding to a 
linear size of 5.2 kpc. The brightest line, [0 11] 3727, has 
a very asymmetric blue wing, as shown in Figure 1111 We 
modelled the line with two gaussians, a narrow component 
(with width Aw = 450kms~^ FWHM) at a velocity very 
similar to the rest-frame velocity, and a broad component 
(All = 830 kms"^ FWHM) blueshifted by 170 kms^V The 
velocity of this component is shown as the large red circle 
in Figure [12] 
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Figure 11. Line profile fitted to the [Oil] and [Olll] lines ex- 
tracted from the central two rows of spectrum A. The dotted 
lines show the two separate Gaussian components which were fit- 
ted. The line separations and strengths of the [O III] doublet were 
constrained to those known from atomic physics, and we required 
the separation of the broad and narrow components to be the 
same as for the [O ll] line. 



The [Olll] 4959, 5007 doublet also showed an asym- 
metric blue wing in the nuclear spectrum, though the broad 
component was noticeably weaker, compared to the narrow 
component, than in the [O ll] 3727 line. We fitted two Gaus- 
sian components to each line, constraining the wavelength 
ratio to the known value and the flux ratios to be 3.0:1, given 
by the transition probabilities; we also constrained the corre- 
sponding components (broad and narrow) to have the same 
line width in each line of the doublet. The widths of the 
gaussians were consistent with the line widths from the fit 
to the [O ll] line, so we constrained them to be the same, al- 
lowing only the line intensities to vary. The resulting fit was 
also acceptable, though with a much weaker broad compo- 
nent. None of the other emission lines — H/?, [Ne ill], [O l] — 
showed any evidence for the presence of a broad component, 
so if present it must be extremely weak. 

To investigate the gas kinematics away from the nu- 
cleus, we extracted spectra from different positions along 
the slit, as described in i; 14.2.31 Except for the nuclear spec- 
trum, there was no evidence for a broad component, so we 
fitted a single gaussian to each line. The resulting velocity 
profiles are shown in Figure [T^ 



4-3.2 lonisation state 

The ionisation state of the MRC B1221-423 is quite low; 
[O ll] 3727 is stronger than [O ill] 5007, while He ll 4686 is ex- 
tremely weak. The line-ratios are characteristic of "low ion- 
isation nuclear emission-line regions" or LINERs (i Heckmanl 
Ll980). The temperature-sensitive [Olll] line intensity ratio 
(A 5007-1- A 4959) /A 4363 is relatively low, indicative of high 
temperatures and/or densities. 

The Ha to H/3 ratio is extreme; we measure Fua/Fup ~ 
6.36, mu ch la rger than the case B recombination value of 
2.86 (e.g. lOsterb r ock i 19891 ). If this were due to reddening, it 
would imply Ay = 2.3 mag, which is much higher than the 
value we derived from the stellar continuum. Ay ~ 0.4 mag 
(§ I4.2.2|) . Interestingly, this is similar to our earlier esti- 
mate of the reddening in the nucleus from modelling the 
spectral energy distributions (Paper I), where we found 
Av = 1.9 mag. If the Balmer emission is mainly from the 
nucleus, we could be seeing the effect of very localised dust. 




-10 1 

Distance from nucleus fkpc) 



Figure 12. Rest-frame velocity profiles along spectrum A, for the 
stellar continuum (solid line) and three emission lines: [Oil] 3727, 
[Olll] 4959, 5007 and [O l] 6300. The positive direction is to the 
south-east. Velocity zero is taken as the velocity of the absorption- 
line component at the position of the nucleus (corrected to the 
heliocentric frame, see ij 14. 2.31 1. The large filled triangle shows the 
velocity of the narrow Hi absorption component (§|3ll, while the 
large filled circle shows the velocity of the broad [O ll] emission. 
The hatched bar shows the full width (FWZI) of the H I absorp- 
tion. 



On the other hand, this ratio could be being distorted by the 
presence of H/3 in absorption; for comparison, we measure 
Fuj/Fup = 0.48, compared to an expected value of 0.47 for 
T — 10'* K. Large ratios of Ha/H/? are seen in many LINER 
sources (|Heckmanlll98d ). 

We can investigate the changing ionisation state of the 
gas by comparing emission-l ine intensity ratios at differ- 
ent positions along the slit (JBaldwin et al.Ml981h . In Fig- 
ure [13] we show the diagnostic diagram for the line ratios 
[O ll]3727/H/3 and [O in]5007/H^, using the gaussian fits to 
the lines described i n i; 14. 3. II We also show the classification 
regions of lLamareille. (2010. '). The diagram clearly shows the 
changing ionisation state at different distances from the nu- 
cleus. Close to the nucleus the line-ratios are consistent with 
ionisation by the AGN, while at larger distances the line ra- 
tios look more like star-forming galaxies. 



5 DISCUSSION AND CONCLUSIONS 

We have presented detailed radio and optical spectroscopic 
observations of MRC B1221— 423, a young CSS source which 
is currently undergoing a merger. Our main findings are 

• a highly distorted radio morphology 

• a large Faraday depth and significant depolarisation in 
the northern, less distorted, lobe of the source 

• the presence of redshifted H l in absorption 

• a range of stellar ages in the galaxy 

• strong emission lines with relatively low ionisation. 

The radio observations of MRC B1221-423 show dra- 
matic interaction between the southern radio jet and the 
host galaxy, together with infalling gas. This jet appears to 
be bent through a full 180°. The alignment of the magnetic 
field along the jet supports the interpretation of the fea- 
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log([0 ll]37Z7/H/5) 

Figure 13. lonisation mechanism diagnostic plot using the Une 
ratios [On]3727/H/3 and [Olll]5007/H/3. Data points represent 
the observed spectra extracted in non-overlapping 2-pixel steps 
along the slit in spectrum A. The red lines show the classification 
regions from Lamareille (2010). The emission lines are AGN-like 
close to the core, and become more star-formation-like further 
out. 



ture as a single connected structure, with the magnetic field 
being oriented circumferentially where the extreme bending 
occurs. This is suggestive of a strong interaction with high 
density material. Interestingly, the RM synthesis shows no 
evidence of a dense interstellar medium where the jet bends. 

The much shorter extent of the radio emission to the 
north, as well as the steep gradient in Faraday depth from 
north to south, suggests the presence of a dense and pos- 
sibly highly magnetised medium local to the source. One 
possibility is that we could be seeing an orientation effect, 
with the northern jet pointed away from us and seen behind 
the galaxy. Large changes of jet direction in radio galaxies 
have been modeled as the result of a helical jet (see e.g. 
iHunstead et al.lll984 [Conwav fc MurphvUiggsl ): it is possi- 
ble that a helix viewed almost head-on could explain the 
extreme change in projected jet direction that we see in 
MRC B1221-423. The lack of evidence for a dense ISM at 
the bend of the jet from the RM synthesis would be consis- 
tent with this interpretation. Alternately, the infalling gas 
could be interacting with the jet, producing both the change 
in jet direction and th e different polarisation in the lobes. 
ISaikia fc Guptal (|2003l ) found a high degree of polarisation 
asymmetry in CSS sources, which they ascribed to interac- 
tion of the jets with infalling material. 

H I was clearly detected in absorption, with a narrow 
component at a velocity +250 kms"'^ compared to the ve- 
locity of the galaxy determined from the stellar component. 
Both this absorption and the radio continuum observations 
suggest the source is in a gas-rich environment. There is also 
a broader wing of H I absorption extending to the systemic 
velocity. The velocity shift with respect to the stellar red- 
shift suggests large amounts of infalling gas. It should be 
noted, however, that the velocities of the Hi, particularly 
of the broad component, are quite similar to the range of 
velocities shown by the emission lines (see Fig. I12p , so it is 



possible the Hi absorption is just part of the overall distri- 
bution of gas. Redshifted Hi absorption is not common i n 
similar sources: of 14 sources studied bv lHolt et al.l l|200a ). 
most showed significantly blueshifted H I absorption, reflect- 
ing outflow. Only two of their sources — PKS B0023— 263 
and PKS B1934-638 — showed redshifted Hi absorption. 

The location of the absorbing gas is not clear. A com- 
parison can be drawn with the radio galaxy NGC 315, which 
shows two absorption systems against the central region 
IjMorganti et al.l [2009|) . In addition to the redshifted nar- 
row component [v ~ -1-490 kms~^), NGC 315 has a broader 
(FWZI ~ 150 kms~^) component re dshifted by ^ 80 km s~^ 
compared to the systemic velocity. iMorganti et al.l (120091 ) 
concluded from velocity gradients seen in VLBI data that 
the broad component arises from gas falling into the nu- 
cleus, while the narrow component arises from more distant 
gas clouds falling towards the galaxy. By analogy, we might 
be seeing a similar situation in MRC B1221— 423. The nar- 
row absorption component could represent a discrete cloud 
at some distance from the galaxy, while the broad absorp- 
tion arises much closer to the nucleus, either from a cir- 
cumnuclear torus, or from the infalling gas that is fueling 
the nuclear activity. Higher resolution Hi observations are 
needed to distinguish these possibilities. 

Optical spectra of MRC B1221-423 reveal disturbed 
kinematics. Most of the ionised gas is at velocities simi- 
lar to that of the stars, but in the nucleus we detect a 
broad blue-shifted component to the emission lines with 
Av = 830kms~^ FWHM, which we interpret as evidence 
for outflow of emission-line gas in the nucleus. Since it is only 
visible in the nucleus, this could arise from gas interacting 
with the radio jet: the fact that it is blue-shifted might sup- 
port the interpretation that the jet axis is strongly inclined 
to the line of sight. 

The ionisation state of the gas is relatively low, with 
[O ll] 3727 stronger than [O ill] 5007, and a very weak 
He II 4686 line. The line ratios change as a function of posi- 
tion, with the emission being dominated by the AGN close 
to the nucleus, and becoming more star-formation-like fur- 
ther out. A similar result was found in t h e sam ple of com- 
pact radio sources studied bv lHolt et al.l (|2009l ) . where the 
ionisation in the nuclear regions is dominated by the cen- 
tral AGN. Like them, we find a lower [O lll] 5007 luminosity 
than extended sources with similar radio luminosity in their 
comparison sample. 

Star-formation activity is evident over the whole sys- 
tem, including the companion galaxy. The youngest stellar 
population is close to the nucleus, consistent with the idea 
that the same inflows of gas that have triggered star for- 
mation are also triggering the activity of the central black 
hole. 

Our previous work on MRC B1221— 423 suggested the 
following scenario: tidal interaction with the companion 
galaxy triggered different phases of star formation in the 
two galaxies, and then, after a substantial delay, the cen- 
tral AGN was triggered to produce the radio emission. The 
current observations support this model. We have found ev- 
idence of past episodes of star formation, probable evidence 
for infalling gas, and indication of interaction between the 
radio jet and its environment. 

Our work confirms the fact that young, powerful ra- 
dio sources tend to live in a gas-rich nuclear environment. 
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Further, we show that the AGN is younger than the star- 
burst. The ages of the stellar populations are comparable 
with the orbital period of the companion, which we esti- 
mate to be ~ 10* yr, based on the projected separation of 
the companion galaxy from the nucleus, together with an 
estimate of the mass of the galaxy fro m the A"-band lumi - 
nosity. The faint arcs to the north-east (jSafouris et al.ll2003l ) 
suggest there have been previous tidal interactions. These 
timescales suggest a time delay oi^ a. few x 100 Myr between 
the tidal event (that triggered star formation) and the on- 
set of AGN activity. S uch delays have been found in several 
nearby AGN (see e.g. J Tadhunter et al.l 1 199a : lEmonts et al.l 
120061 : iDavies et al.ll2007l ). but this is possibly one of the best 



examples. In recent work, IWild et al.l (|20ld ) constructed a 
sample of starburst galaxies from the Sloan Digital Sky Sur- 
vey. They found that the rate of accretion onto the central 
black hole, as measured by the strength of the [O ill] line, 
rises steeply about 250 Myr after the onset of star forma- 
tion. They speculate that the accretion onto the black hole 
may be being regulated by the starburst, with the accretion 
being dominated by low-velocity stellar winds, and accretion 
at early times being suppressed by supernovae. 

We have also found indication of a possible infall of cold 
gas that could provide the fuel for the AGN. We note that 
the age distribution of the stellar populations is consistent 
with a picture in which cold or cooling gas is driven by the 
merger towards the centre of the host galaxy to fuel the AGN 
from relatively large radii, and that this gas begins to form 
stars while it is still well away from the nucleus. This con- 
trasts with the situation in lower-luminosity radio sources in 
cluster environments, which, even though they sometimes 
contain the most massive black holes, are accreting in an 
environment dominated by hot gas and are associated with 
an accretion mode that produces low-power radio jets. 
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APPENDIX A: FLUX MEASUREMENTS 

In Table lAll we present radio flux measurements of 
B1221--423 from the literature. 
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Table Al. B1221-423 radio flux densities. 



Frequency 
(MHz) 



(mJy) 



Survey/ 
Telescope 



Reference 



80 


8000 ± 1300° 


CCA 


160 


8100 ± 800" 


CCA 


408 


5080 ± 130 


MRC 


468 


4670 ± 810 


PKS 


635 


3910 ± 270 


PKS 


843 


3300 ± 200 


MOST 


843 


3110 ±90 


SUMSS 


960 


3070 ± 170 


PKS 


1344 


2360 ± 50 


ATCA 


1382 


2290 ± 70 


ATCA 


1410 


2470 ± 130 


PKS 


2280 


1650 ± 80 


ATCA 


2496 


1600 ± 50 


ATCA 


2650 


1700 ± 50 


PKS 


2700 


1620 ± 50 


PKS 


4740 


1050 ± 50'' 


ATCA 


4790 


1010 ± 50'"= 


ATCA 


4800 


1010 ± 20 


ATCA 


4800 


999 ± 41 


ATCA 


4850 


1090 ± 60"* 


PMN 


4850 


1140 ±60 


PMNM 


5009 


1000 ± 50 


PKS 


8400 


630 ± 70 


PKS 


8640 


610 ± 20 


ATCA 


18 500 


357 ±18 


ATCA 


B 880.5 


310 ± 30 


ATCA 


3 648.5 


290 ± 30 


ATCA 


19 904 


330 ± 17 


AT20G 


22 000 


347 ± 35 


ATCA 



Slee 



Slee fc HiKKins ( 1973) 

Slee (1977, 1995) 
Large ct al. (1981) 
Wills (1975) 
Wills (1975) 

^urgess_&_jiunstead (2006a) 
Mauch et aL (2003) : Murphy et al. 



Wills (1975J 
This paper 
This paper 
Wills (1975) 
This paper 
This paper 
Wills (1975) 

Bolton fc Shimmins (1973) 
airgess & Hunstead (2006b) 
D uncan fc Sproats (1992 ) 
Safouris et al. (2003) 
^ur£h^^t_al. 
^right_^t_^l^ 
Gregory et al. 
Wills (1975) 

Wright et aL (1990) : Wright fc Otrupcek 
Safouris et aL (2003) 
Ricci et al. (2006) 
This paper 
This paper 
Murph y et al. (2010) 
Ricci et al. (2006) 




Notes: CCA - Culgoora Circular Array; MRC - Molonglo Reference Catalogue; MOST - Mo- 
longlo Observatory Synthesis Telescope; SUMSS - Sydney University Molonglo Sky Survey; 
PKS - Parkes Catalogue; ATCA - Australia Telescope Compact Array; PMN - Parkes- 
MIT-NRAO Survey; PMNM - PMN Map Catalogue; AT20G - Australia Telescope 20 GHz 
Survey 

''Flux density uncertainty calculated using the error analysis from the original CCA cata- 
logue, but not inclu ding the original uncer tainty associated with the flux density scale; the 
scale was refined bv lSlee fc SiegmanI l ll98Sr . 

''Flux density uncertainty not stated in the reference; a ±5 per cent uncertainty has been 
assumed. 

'^The original measurement has been rescaled so that it is consistent with the current ATCA 
flux density scale. 
** General- width fit; the flxed-width fit has a similar fiux density (1100 ± 60 mJy). 
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